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Summary 

A 22-percent dynamically scaled replica of the 
X-29A forward-swept- wing airplane has been flown in 
radio-controlled drop tests at the Langley Plumtree 
Test Site. Flight data, which were recorded from 
early flights in the test program, consisted mainly of 
large-amplitude maneuvers over wide angle-of-attack 
ranges with several uncontrolled wing rock episodes. 
A system identification study of the recorded data 
was undertaken to examine the stability and control 
derivatives that influence the lateral behavior of this 
vehicle with particular emphasis on an observed wing 
rock phenomenon. All major lateral stability deriva- 
tives and the damping-in-roll derivative were identi- 
fied for angles of attack from 5° to 80° by using a 
data-partitioning methodology and a modified step- 
wise regression algorithm. No control effectiveness 
derivatives could be identified from these first flights 
in the test program. 

Introduction 

Tests of dynamically scaled airplane models con- 
tinue to be the most reliable source of information on 
high-angle-of-attack, flight dynamic characteristics 
prior to the actual flight tests of a particular airplane. 
A properly scaled dynamic model may be thought of 
as a simulator with the proper values of the various 
aerodynamic and inertial parameters. One such dy- 
namic model test method, the radio-controlled drop 
model technique, involves dropping an unpowered, 
dynamically scaled model from a helicopter, flying it 
through a series of test maneuvers by remote con- 
trol from the ground, and finally recovering it with 
a parachute. These tests can provide an important 
step toward understanding the dynamic characteris- 
tics of an airplane configuration since the testing is 
done in the actual flight environment. In addition, 
these models do not endanger a human pilot in an ex- 
pensive airplane prototype. The emergence of good, 
miniature, instrumentation packages and telemetry 
systems makes it possible to measure and transmit 
air data, control surface deflections, and airplane re- 
sponses to a ground station where they are recorded 
for later analysis. 

Drop tests of a model of the X-29A forward- 
swept-wing airplane were conducted by the Langley 
Low-Speed Aerodynamics Division (LSAD). A pur- 
pose of these tests was to verify a wind-tunnel pre- 
diction of wing rock for angles of attack greater than 
20°. In the first two data flights, violent lateral os- 
cillations occurred in this angle-of-attack range, thus 
culminating in a roll departure. In order to ascer- 
tain the source of these severe rolling motions, the 
recorded data from those tests were provided to the 


Langley Guidance and Control Division where a sys- 
tem identification analysis was conducted. The pur- 
pose of this paper is to present the results of that 
analysis. A thorough description of the drop tests is 
presented in references 1 and 2. 

First, the conduct of the drop test procedure and 
the system identification methodology that was ap- 
plied to the flight data are described. Results of the 
analysis in the form of lateral stability and damping 
derivative estimates are then discussed and compared 
with wind-tunnel predictions from reference 3. Re- 
sults of a validation scheme for checking the integrity 
of the estimates are also presented and discussed. A 
concluding section summarizes the most important 
results. 

Symbols and Abbreviations 

ay lateral acceleration, g units 

b wing span 

C a general aerodynamic force and 

moment coefficient 

Ci rolling-moment coefficient, Mx/qSb 

C n yawing-moment coefficient, Mz/qSb 

Cy lateral-force coefficient, Fy/qS 

c wing mean aerodynamic chord 

Fy force along lateral body axis 

g acceleration due to gravity 

(lg « 32.174 ft/sec 2 ) 

lx, ly, Iz moments of inertia about longitu- 

dinal, lateral, and vertical body 
axes 

J cost function 

t + 1 number of parameters in regression 

equation 

Mjf, Mz rolling and yawing moments 

N number of data points in bin 

p roll rate, rad/sec or deg/sec 

q dynamic pressure, pV 2 / 2 

S wing area 

t time, sec 

V free-stream velocity 

X{ independent variable (regressor) in 

regression equation 

y dependent variable (aerodynamic 

force or moment coefficient) in 
regression equation 


a angle of attack, deg or rad 

d rate of change of angle of attack, 

deg/sec or rad/sec 

/ 3 sideslip angle, deg or rad 

6 a aileron deflection angle, deg or rad 

6 r rudder deflection angle, deg or rad 

6i nondimensional stability or control 

derivative in regression equation 

a standard deviation 


Subscript: 

0 bias term in regression equation 

Abbreviations: 


LSAD 

MSR 

SAS 

TED 

TEU 


Low-Speed Aerodynamics Division 
modified stepwise regression 
stability augmentation system 
trailing edge down 
trailing edge up 


A dot ( ) over a symbol denotes a derivative 
with respect to time. A circumflex ( A ) denotes an 
estimated value. 

Derivatives of aerodynamic coefficients C a (a = 
Y, /, n) referenced to a system of body axes with the 
origin at the airplane center of gravity: 
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Flight Vehicle Description and Drop Test 
Procedure 

Flight data were generated during drop tests of 
an unpowered, 22-percent dynamically scaled replica 
of the forward-swept-wing X-29A research airplane. 
The model is equipped with controllable flaperons, 
rudder, canard, and strake flaps as depicted in fig- 
ure 1. The flaperons act both symmetrically and 
differentially for lift and roll, respectively. Canard 
and strake flap deflections are purely symmetric for 
pitch control. The model is unpowered with the re- 
sult that no fuel is burned and the center-of-gravity 
location is constant throughout the flight. A sum- 
mary of the mass and geometry characteristics of the 
model is presented in table I. 


Drop tests are conducted at the Langley Plumtree 
Test Site. A schematic representation of the drop 
test procedure is presented in figure 2. (See ref. 1.) 
A helicopter carries the model to an altitude of ap- 
proximately 6000 ft and then releases it for a radio- 
controlled descent. The overall model control loop 
consists of downlinked flight data, a digital com- 
puter on which the control laws are implemented, 
and a command uplink. The pilot, who sits at a 
control station monitoring data displays and a video 
image of the model, provides three-axis commands 
through a customized control box. The instrumen- 
tation package of the model measures angle of at- 
tack, angle of sideslip, airspeed, linear accelerations, 
angular rates, and control surface deflections. Data 
time histories are recorded at the rate of 100 sam- 
ples per second. Recorded angular rates are numer- 
ically differentiated after the flight to obtain angu- 
lar accelerations. Flights are numbered sequentially 
by LSAD. This report presents parameter estimates 
based on data from flights 2, 3, 5, 6, 7, and 8. Data 
from flight 4 were obtained, but that flight consisted 
mainly of a well-developed spin; thus, the informa- 
tion content of those data was insufficient for the es- 
timation of lateral stability parameters. 

Flights 2 and 3 operated with a longitudinal 
stability augmentation system (SAS) in effect and 
with virtually no lateral augmentation. These flights 
covered an angle-of-attack range from 5° to about 
35°. With essentially no lateral SAS, a sudden onset 
of poorly damped, uncommanded rolling oscillations 
(wing rock) occurred as the model entered certain 
angle-of-attack ranges. The wing rock behavior of 
the model, which made control difficult at angles 
of attack greater than 20°, eventually resulted in a 
roll departure. The first set of parameter estimates 
presented in this paper is based on data from flights 2 
and 3 and covers angles of attack up to 32°. 

Flights 5, 6, 7, and 8 were conducted with longi- 
tudinal and lateral SAS operating, which permitted 
control of the vehicle to much higher angles of at- 
tack. Wing rock oscillations were suppressed, and 
data were produced at angles of attack up to 80° 
during pull-up maneuvers. The second set of param- 
eter estimates presented in this paper is based on 
data from these flights for angles of attack between 
30° and 80°. None of the flight data analyzed in 
this report consisted of maneuvers designed specifi- 
cally for parameter estimation purposes. Therefore, 
the flights did not contain simple control doublets 
and pulses that are necessary for a straightforward 
estimation of control effectiveness derivatives. No 
control effectiveness derivatives could be consistently 
identified from these flights. A summary of several 
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important characteristics associated with each flight 
is presented in table II. 

Analytical Methods 

A system identification study typically addresses 
two problems. The first problem involves model 
structure determination, and the second problem in- 
volves estimation of the parameters associated with 
terms included in the model. The system identifica- 
tion study described in this report utilized a regres- 
sion technique to obtain parameter estimates from 
flight data. Prior to the regression procedure, a par- 
titioning process was applied to the flight data. 

Data Partitioning 

The X-29A drop tests contained many transient 
and large-amplitude maneuvers over a wide angle-of- 
attack range. In order to produce stability parameter 
estimates that reflect the overall variation of the aero- 
dynamic coefficients with angle of attack (a), a data- 
partitioning technique was employed. This technique 
permits the analysis of unplanned or uncontrolled 
large-amplitude maneuvers such as those experienced 
in flights 2 and 3. Data partitioning requires the 
sorting of data points into angle-of-attack bins, each 
of which covers a small angle-of-attack range (usu- 
ally only 1° or 2°). A single bin will contain data 
from all portions of the maneuver during which the 
model was in that particular angle-of-attack range. 
These portions of flight data do not have to be, and 
often are not, contiguous in time. Figure 3 illus- 
trates the manner in which a single maneuver may be 
partitioned into several bins on angle-of-attack space 
(a space). This process is thoroughly described and 
demonstrated in reference 4. 

The regressional analysis (to be discussed below) 
employs a simple curve-fitting algorithm and does 
not require data that are contiguous in time. Since 
continuity in time for data in a particular bin is not 
needed, it is possible to combine data from separate 
maneuvers or even from separate flights into a single 
bin. In this way, information from several flights can 
be analyzed in one data set. The a dependence of 
the data in each bin is removed since the a range 
of each bin is small. Any bin that appears to retain 
some dependence may be broken into smaller bins 
if an adequate number of data points remain in each 
bin. This partitioning procedure, which has generally 
been applied on a space, may also be applied with 
respect to any other variable desired. 

The partitioning of a maneuver that covers a large 
angle-of-attack range in a short period of time tends 
to yield a large number of bins that contain few data 
points. This was particularly true of maneuvers in 


flights 5, 6, 7, and 8 during which the lateral SAS 
was in effect, thus permitting controlled flight over 
a much broader a range than in flights 2 and 3. In 
order to produce a sufficient number of data points 
for regressional analysis on each bin, it was necessary 
to combine like bins from separate maneuvers and 
flights. An example of this is shown in figure 4 
which illustrates the manner in which two separate 
maneuvers from flights 2 and 3 produced a data 
ensemble that contained sufficient points to create 2° 
bin sizes on an angle-of-attack range from 19° to 31°. 
Maneuvers from flights 5, 6, 7, and 8 were combined 
to produce an ensemble of 3° bin sizes on an a range 
from 30° to 80°. Estimates for lateral-directional 
stability and roll damping on these two a ranges will 
be presented and discussed in this report. 

Modified Stepwise Regression 

Once the partitioning has been completed, the pa- 
rameter estimation may be conducted for each bin. 
In this study, an equation error method referred to 
as modified stepwise regression (MSR) was employed 
as a means of model structure determination and pa- 
rameter estimation. The MSR algorithm (ref. 5) is a 
modification of linear regression which can approxi- 
mate the structure of aerodynamic model equations 
and estimate the model parameters. When the analy- 
sis is conducted separately for each bin of partitioned 
data, the MSR yields a series of parameter estimates 
that reveal trends with respect to angle of attack. 
MSR successively enters a new independent variable 
(regressor) one at a time into a regression equation. 
Each new regressor is selected from a list of candidate 
variables established by the user before initiation of 
the regression procedure. 

The general form of a regression equation for an 
aerodynamic model can be expressed as the poly- 
nomial expansion 

y(t) = #o + #1 ^1 W + . . . + 0n z n (£) 

where y(t) represents an aerodynamic force or mo- 
ment coefficient, xj(t) represents a regressor consist- 
ing of an airplane output or control variable (or some 
combination thereof), and 0j represents the stabil- 
ity or control derivative associated with the regres- 
sor Xj(t). The new regressor is selected from the list 
of candidate variables based on the correlation of its 
time history with that of the aerodynamic coefficient 
y(t). Values for parameters 6j are obtained by mini- 
mizing the cost function J$r defined as 

n r e l 2 

Jsr = e j *>(*) 

i=i _ j= i 
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where N equals the number of data points in the 
bin being analyzed and (£ + 1) is the number of pa- 
rameters in the current regression equation. This 
cost function represents the least-squares difference 
between the measured value of the aerodynamic co- 
efficient (based on sensor data) and the predicted 
value of the aerodynamic coefficient (based on the 
current regression model). The “best” parameter es- 
timates for 9j will be those that minimize this cost 
function. The iterative process of adding regres- 
sors to the aerodynamic model, and then estimat- 
ing the new parameters associated with those regres- 
sors (and the ones previously in the model), allows 
one to assess the impact of each new addition to the 
model and its ability to predict the given aerody- 
namic coefficient. In each step of the process, all 
terms in the regression are evaluated for their sig- 
nificance in minimizing the cost function, and those 
found insignificant according to a statistical selec- 
tion criterion are deleted. The process continues un- 
til no new regressors are selected from the list of 
candidate variables. In applications to actual flight 
data, the final model produced in this manner may 
be overparameterized; thus, criteria for selecting an 
adequate model from an earlier step in the regres- 
sion have been developed and are discussed in refer- 
ence 6. A more detailed description of modified step- 
wise regression may be found in references 5, 6, and 
7, which also include information relating to the use 
of this algorithm on a digital computer. 

Results and Discussion 

Estimates of lateral-directional stability deriva- 
tives for the X-29A drop model have been obtained 
from the six-flight data ensemble that was parti- 
tioned on a space. The combined estimates cover 
an a range from 5° to 80°. Estimates from flights 2 
and 3, which cover an a range from 5° to 35°, are 
presented in figures 5 and 6 (solid symbols). Open 
symbols, for comparison, represent wind-tunnel data 
from a 16-percent-scale model of the X-29A that 
was tested in the Langley 30- by 60-Foot Tunnel 
(ref. 3). The wind-tunnel data used for comparison 
were taken for the model with control settings of the 
canard at -25°, the flap at 17.5°, and the strake at 
30°. Several discrepancies are seen between wind- 
tunnel and flight test results for a > 20°. These 
may arise because the dynamics of the flight maneu- 
vers give rise to time-dependent phenomena that are 
not evident in static wind-tunnel tests. Evidence of 
similar discrepancies (although both are correct mea- 
surements for their respective tests) is given in ref- 
erences 8 and 9. The estimates for C n(3 confirm the 
wind-tunnel prediction that directional stability is re- 
tained to at least a = 30°. Since wing rock episodes 


occurred at angles of attack as low as 17°, this re- 
sult indicates that loss of directional stability is not 
a causal factor in the initiation of wing rock. A de- 
crease in directional stability is observed for a = 15° 
and above in figure 5. This drop in C n 0 may arise 
from a blanketing effect of the fuselage and wing wake 
upon the vertical fin in this angle-of-attack region. 
The estimates for show a slightly lower (but still 
stable) negative dihedral effect than was predicted by 
the wind tunnel. 

Estimates for roll damping in figure 6 indicate a 
region of unfavorable damping for angles of attack 
between 17° and 24°. Another region of unfavor- 
able damping exists for angles of attack greater than 
27°. The existence of unfavorable roll damping is a 
likely cause of the wing rock oscillations observed in 
these a ranges. Wind-tunnel investigations have re- 
vealed that the loss of roll damping is caused by a 
pronounced asymmetry of fuselage forebody vortices 
under a rolling condition. (See ref. 3.) It is interest- 
ing to note that C ^ remains negative throughout the 
regions of unfavorable roll damping, thus indicating 
that wing rock is a dynamic instability manifested by 
a statically stable system. 

A close examination of the data from flights 2 
and 3 reveals that two angle-of-attack regions of dif- 
fering wing rock behavior can be defined. By plot- 
ting sideslip angle (3 versus roll rate p for wing rock 
episodes from these two flights, it becomes apparent 
that the oscillations are limit cycles at certain times 
and divergent cycles at others. The limit-cycle oscil- 
lations are indicated by the closed circular patterns 
at the center of the phase plane plots in figure 7. 
Divergent oscillations are indicated by the outward 
expanding spirals. Each segment of wing rock data 
shown in figure 7 exhibits these two regions of oscilla- 
tion. Examination of these patterns revealed that all 
limit-cycle oscillations occurred at angles of attack 
between 17° and 27°. Above 27°, the wing rock os- 
cillations are divergent. The roll-damping estimates 
from figure 6 reinforce the existence of these trends. 
Because essentially no lateral SAS was operational 
in flights 2 and 3, these data represent the “natural” 
behavior of the vehicle with no control action taken 
to suppress the wing rock oscillation. 

Estimates of lateral-directional stability deriva- 
tives for a — 30° to 80° are presented in figures 8 
and 9. These estimates (solid symbols) are based 
on flights 5, 6, 7, and 8. A constant offset of 
approximately 0.6 is observed between wind-tunnel 
and flight estimates for Cy r The character of the 
two curves is identical, although the constant offset 
persists throughout the angle-of-attack range. The 
development of this offset may be observed in the 
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estimates of Cy g for a = 5° to 35° in figure 5. The 
offset starts at a « 24° and reaches a value of about 
0.6 at a = 32°. It then persists throughout the re- 
mainder of angles of attack shown in figure 8. The 
appearances of this offset in two independent sets of 
parameter estimates tend to corroborate one another. 
No cause was found for this discrepancy. 

The estimates for C[ 0 on a = 30° to 80° show 
some offset from wind-tunnel predictions. Negative 
(favorable) dihedral effect is confirmed up to a = 80°. 
The flight data indicate a slightly less negative value 
than was observed in the wind tunnel for a = 30° to 
65°, but the data agree well for a = 65° to 80°. 

The estimates for C nf) at a — 30° to 80° based on 
the drop tests indicate a loss of directional stability at 
a = 34° (fig. 8). This may represent the culmination 
of the blanketing effect on the vertical fin that was 
observed to begin at about a = 14° in figure 5. The 
flight estimates indicate less directional stability than 
was observed in the wind tunnel. Wind-tunnel data 
indicate a reestablishment of static directional stabil- 
ity above a = 40°, whereas flight estimates indicate 
that directional instability persists; however, a slight 
trend toward recovery is observed between a = 50° 
and 65°. Reestablishment of directional stability in 
the wind-tunnel tests was credited to the unusual 
vortex aerodynamics associated with the long slender 
forebody and flat elliptical cross section of the X-29A 
at high angle of attack. (See ref. 3.) This factor may 
be the cause of the trend toward C nf) recovery be- 
tween a = 50° and 65° that was observed in the 
estimates of drop test parameters. If so, the effect 
does not appear to be as pronounced in flight as 
was observed in the wind tunnel. It is possible that 
the drop model does not maintain static conditions 
long enough for the fuselage vortices to develop com- 
pletely and take their full effect. This discrepancy 
between wind tunnel and flight estimates for direc- 
tional stability in the poststall regime should be in- 
vestigated in subsequent flights of the X-29A drop 
model. 

Since many of the maneuvers in flights 5, 6, 7, 
and 8 consisted of rapid pull-ups through the regions 
of unfavorable roll damping, the dynamics associated 
with changing angle of attack were considered in pro- 
ducing the parameter estimates. To include a consid- 
eration of these dynamics, the flight data were par- 
titioned with respect to the sign of a as well as the 
magnitude of a. The more selective bins contained 
fewer data points, and therefore they contained 
less information for the regression using MSR. Sub- 
sequent application of the modified stepwise regres- 
sion produced a slightly better fit of the flight data 
than did previous estimates. Estimates of C/ p for 


a > 0 and d < 0 from flight data are compared with 
wind-tunnel predictions of C\ in figure 9. The values 
of C{ p for a > 0 indicate less favorable roll damping 
than for d < 0. Wind-tunnel predictions lie roughly 
between the two sets of estimates. This disparity be- 
tween estimates of C/ for d > 0 and d < 0 in the 
angle-of-attack range between 30° and 80° indicates 
that dynamic conditions are giving rise to hysteresis 
in the vehicle behavior with respect to angle of at- 
tack. The estimates of a less favorable (7/ for d > 0 
than for d < 0 suggest that the vehicle is more prone 
to wing rock episodes while a is increasing through 
the critical angle-of-attack ranges rather than de- 
creasing. This observation of hysteresis with d is of 
particular significance, since it indicates that the pa- 
rameter estimation techniques employed herein may 
be capable of quantifying the effects of a somewhat 
elusive phenomenon such as dynamic stall if given 
adequate flight test data. A greater tendency toward 
wing rock was noted in the drop tests during pitch-up 
segments of maneuvers in the higher angle-of-attack 
ranges. Estimates of the static lateral-directional sta- 
bility derivatives Cy 0 , C{ 0 , and C H0 were essentially 
identical for the conditions of d > 0 and d < 0. 
Tables III to V present a numerical listing of the es- 
timates of stability derivatives from figures 5 to 9. A 
summary of the estimates of lateral stability param- 
eters for the angle-of-attack range from 0° to 80° is 
presented in figures 10 and 11. 

Validation 

A good check on the validity of parameter esti- 
mates obtained from a system identification study is 
to use the estimates to predict the forces and mo- 
ments which act on the flight vehicle. The following 
section describes the manner in which predicted time 
histories of the drop model behavior based on the 
parameter estimates were compared with predictions 
based on wind-tunnel data. Both sets of predictions 
were plotted against actual time histories of the ve- 
hicle behavior based on telemetered flight data. 

Method 

Time histories for the coefficients Cy, Ci, and 
C n have been computed and compared as a means of 
checking the integrity of parameter estimates. Three 
sets of time histories were computed for each coef- 
ficient. First, time histories of force and moment 
coefficients were computed from flight data using the 
telemetered measurements from accelerometers and 
differentiated rate gyros onboard the vehicle. These 
values represent measurements of the actual forces 
and moments experienced by the model during a 
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flight. Time histories for the force and moment coef- 
ficients were also computed based on simplified aero- 
dynamic models using wind-tunnel estimates of the 
lateral stability parameters. Coefficient time histo- 
ries were then computed from the same simplified 
aerodynamic models using drop test estimates of the 
lateral stability parameters that were generated by 
the MSR algorithm with data partitioning. The sim- 
plified models used to compute the side-force and mo- 
ment coefficients are as follows: 

C Y {t) = C Yt o + C Y0 [a(t)\ m + C Y6r [a(*)]«r(t) 

C t (t) = C it o + C l0 [a(t)} m + C lp [a(t))p(t)± 

+ Ci 6 a [a{t)]S a (t) 

Cn(t) = Cn,o + Cfi/jfaW] /^(0 "I" [ Q! (^)]^r(^) 

Note that the stability coefficients are functions of 
a(t) and are obtained from interpolation of the pa- 
rameter estimate summaries shown in figures 10 
and 11. The control derivative estimates Cy^, Cn Sr , 
and Cu are taken from reference 10, and the same 
estimates are used in calculating both drop test and 
wind-tunnel coefficients. A comparison of time his- 
tories obtained from stability parameter estimates 
with the actual time histories provides a general in- 
dication of the integrity of the parameter estimates. 
An example of a side-force-coefficient time history 
is shown in figure 12. Figure 12(a) compares coef- 
ficient predictions based on wind-tunnel parameter 
estimates with the measured coefficient values (solid 
line). Figure 12(b) compares coefficient predictions 
based on MSR parameter estimates from drop tests 
with the measured coefficient values. Figure 12(c) 
shows the angle-of-attack time histories which corre- 
spond to the coefficient comparisons in figures 12(a) 
and 12(b). These allow the reader to determine the 
original parameter estimates from figures 10 and 11 
that were used to compute the coefficients. 

Figures 12 and 13 present time history compar- 
isons for the side-force coefficient Cy from two sepa- 
rate flight segments. The wind-tunnel estimates ap- 
pear to underpredict the amplitude of the side-force 
oscillations that were observed during these portions 
of the drop tests. Note that the time histories cover 
an angle-of-attack range of approximately 0° to 60° 
in these two particular flight segments. The time 
histories shown indicate that the estimates of system 
identification parameter, which were derived from a 
partitioned data ensemble from several flights, pro- 
vide a good ability to predict the vehicle behavior 
during a particular flight maneuver. 


Figure 14 presents another side-force-coefficient 
time history. This maneuver includes a higher angle- 
of-attack range than the previous two. At approxi- 
mately 11 sec into this maneuver, the angle of attack 
exceeded 50°. Wind-tunnel predictions of the coef- 
ficient Cy apparently produced the wrong sign on 
the side-force oscillations observed in this time seg- 
ment. Reviewing the wind-tunnel estimates of Cy Q 
in figure 10 reveals that a change in the sign of this 
parameter was, in fact, predicted at approximately 
a = 50°. This wind-tunnel prediction of positive 
Cyp apparently does not coincide with the actual be- 
havior of the drop model. The overall agreement be- 
tween the system identification predictions and the 
measured coefficient time history is somewhat poorer 
than in the two previous maneuvers, thus suggesting 
that the simplified model being used for Cy is less 
adequate at these higher angles of attack. 

Figures 15 and 16 present time histories for 
the rolling- moment coefficient C/. It appears that 
the system identification predictions provide better 
agreement with the measured coefficient values at the 
peaks of the rolling-moment oscillations. This im- 
provement is particularly apparent in figure 16 during 
the interval between 13 and 18 sec. Angles of attack 
during this interval lie within the range wherein data 
were partitioned with respect to the sign of a as well 
as to the magnitude of a. This distinction as to the 
dynamic condition of angle of attack may be respon- 
sible for the improved prediction of rolling moment 
at the peaks of the oscillations. The simplified model 
used appears to provide a relatively good agreeement 
even in the higher angle-of-attack ranges. 

A comparison of yawing-moment time histories is 
presented in figure 17. Most of the flights contained 
poor yaw excitation. The primary parameter influ- 
encing the yawing coefficient C n is the directional 
stability derivative Cnr For the particular flight seg- 
ment shown in figure 17, it appears that the wind- 
tunnel estimates mispredict the sign of the yawing- 
moment oscillations. The angle-of-attack range for 
which this occurs (approximately 55° to 80°) cor- 
responds with the range in which wind-tunnel esti- 
mates predicted a recovery of directional stability, 
whereas the drop test estimates predicted a contin- 
ued loss of directional stability (fig. 10). Apparently, 
this recovery of directional stability did not occur 
during the flight of the drop model. 

Conclusions 

Estimates of lateral-directional stability param- 
eters and roll-damping parameter have been calcu- 
lated using modified stepwise regression on flight 
data from drop tests of a 22-percent-scaled replica 
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of the X-29A forward-swept-wing research airplane. 
Data partitioning was used to obtain trends in aero- 
dynamic parameters with angle of attack. Data from 
the first two flights were used to obtain parameter 
estimates that covered a range of angle of attack a 
from 5° to 35° and showed relatively good agree- 
ment with wind-tunnel measurements. Retention of 
directional stability was confirmed up to a = 30°, 
and a slightly less negative dihedral effect was ob- 
served than was predicted from wind-tunnel tests. 
Estimates of the variation of rolling-moment coeffi- 
cient with roll rate ( Cj ) indicate a loss of favorable 
roll damping at a = 17° that is followed by a brief 
recovery of damping for a — 24° to 27°; then, un- 
favorable damping occurs again for a > 27°. These 
successive changes in roll-damping character appear 
to give rise to two differing types of wing rock be- 
havior. The first is a limit-cycle oscillation for angles 
of attack between 17° and 27°. The second is a di- 
vergent oscillation at angles of attack greater than 
27°. 

Estimates were also made based on data from four 
flights covering an angle-of-attack range from 30° to 
80°. A constant negative offset from wind-tunnel pre- 
dictions was observed in estimates for the variation of 
side- force coefficient with sideslip angle (Cy^). Esti- 
mates for the variation of rolling- moment coefficient 
with sideslip angle (C/^) showed very good agree- 
ment with wind-tunnel measurements. The config- 
uration is statically stable in roll up to a = 80°. 
Significantly lower directional stability was observed 
than that predicted by wind-tunnel tests. The con- 
figuration loses directional stability at a = 34° and 
remains unstable or marginally stable for all angles 
of attack above this value. Estimates for Ci p were 
produced by partitioning with respect to the sign of 
the rate of change of angle of attack (a) as well as 
to the magnitude of angle of attack. Resulting val- 
ues for roll damping indicate less favorable C\ p for 
a > 0 than for a < 0. Therefore, a greater tendency 
toward wing rock is expected when a is increasing 
rather than decreasing. The dynamic condition of 
angle of attack appears to give rise to this hysteresis 
in the behavior of the X-29A drop model. A vali- 
dation scheme was devised to investigate the ability 
of these parameter estimates to predict the actual 
behavior of the drop model. A better predictive ca- 
pability was found when compared with predictions 
based on wind-tunnel estimates. This validation ex- 
ercise presents the following general conclusions: 

1. Both wind-tunnel measurements and system 
identification estimates provide reasonable predic- 
tions of the forces and moments which act on the 
model in the low-to-mid angle-of-attack range. 


2. Simplified math models used in the validation 
exercise to predict the forces and moments acting on 
the model become less adequate as angle of attack 
becomes larger. 

3. System identification parameter estimates ap- 
pear to reflect the forces and moments which act on 
the drop model at high angles of attack more closely 
than wind-tunnel measurements. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
March 13, 1989 
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Table I. Summary of Mass and Geometric Characteristics of 
X-29A Drop Model 

Weight, lb 276.5 

Moment of inertia: 

I x , slug-ft 2 3.83 

I Y , slug-ft 2 43.62 

Iz, slug-ft 2 43.56 

Wing: 

Span, ft 5.98 

Mean aerodynamic chord, ft 1.58 

Area, ft 2 8.96 

Aspect ratio 4.0 

Sweep (leading edge), deg —29.27 

Sweep (quarter-chord), deg —33.73 

Taper ratio 0.404 

Dihedral angle, deg 0 

Flaperon area (total), ft 2 1.39 

Flaperon deflection, deg 17.5 (TED) to -27.7 (TEU) 

Canard: 

Span, ft 3.0 

Area (total), ft 2 1.79 

Aspect ratio 1.47 

Taper ratio 0.318 

Deflection, deg 30 (TED) to —60 (TEU) 

Vertical tail: 

Height, ft 1.47 

Area, ft 2 1.63 

Aspect ratio 1.33 

Taper ratio 0.306 

Rudder area, ft 2 0.335 

Rudder deflection, deg ±30 

Strake flap: 

Area (total), ft 2 0.504 

Deflection, deg ±30 


Table II. Summary of Flight Characteristics 


Flight 

Range of a , deg, 
on which derivatives 
were identified 

Stability augmentation in effect 

Description of flight 

1 

Not used 

Longitudinal; very weak lateral 

General checkout flight 

2 

5-32 

Longitudinal; very weak lateral 

Shallow pull-ups, turns, 
wing-rock-induced departure 

3 

5-32 

Longitudinal; very weak lateral 

Shallow pull-ups, turns, 
wing-rock-induced departure 

4 

Not used 

Longitudinal; wing rock suppression 

Well-developed spin, 
wing rock suppression 

5 

6 
7 

30-80 

Longitudinal; wing rock suppression; 
spin suppression 

Steep pull-ups, turns, spin 
prevention 
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Table III. Estimates of Modified Stepwise Regression for Lateral-Directional 
Stability Derivatives With Standard Deviations 


Q, 

deg 

C Y 0 

Cl 

0 

Cn 0 

Estimate 

a value 

Estimate 

a value 

Estimate 

cr value 

8 

-0.954 


0.030 

0.002 

0.233 

0.004 

12 

-.848 

mtmm 

.036 

.002 

.266 

.003 

16 

-.807 


-.084 

.003 

.221 

.004 

18 

-.831 


-.073 

.004 

.188 

.003 

20 

-1.050 

wKSmm 

-.102 

.003 

.168 

.003 

22 

-1.350 

.049 

-.099 

.004 

.150 

.004 

24 

—1.170 

0.042 

-0.118 

0.004 

0.113 

0.004 

25 

-.612 

.060 

-.134 

.004 

.088 

.006 

26 

-.700 

.017 

-.142 

.003 

.089 

.004 

27 

-.700 

.020 

-.114 

.003 

.071 

.005 

28 

-.696 

.039 

-.147 

.003 

.073 

.011 

30 

-.720 

.048 

-.134 

.004 

.042 

.016 

32 

-0.740 

0.046 

-0.126 




34 

-.798 

.037 

-.122 



dESS 

37 

-1.060 

.049 

-.153 




40 

-1.110 

.048 

-.120 



.013 

43 

-.989 

.027 

-.145 


-.145 

.012 

46 

-1.070 

.022 

-.210 



.011 

49 

-1.030 


-0.170 




52 

-.560 


-.124 


-.114 


55 

-.365 

.042 

-.127 



.009 

58 

-.110 

.011 

-.114 



.017 

61 

-.273 

.013 

-.116 


1 

.007 

64 

-.144 

.020 

-.132 


1 

.008 

67 

-0.314 


-0.132 



■1H 

70 

-.550 


-.145 




73 

-.584 


-.136 




76 

-.476 

' 1 1 

-.154 




79 

-.267 


-.176 
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Table IV. MSR Estimates of Roll Damping Derivative With Standard 
Deviations — Data Partitioned by Positive and Negative & 


a, 

deg 

Ci p for d > 0 

Ci p for a < 0 

Estimate 

a value 

Estimate 

a value 

8 

-0.491 




12 

-.362 


-.362 


16 

-.191 


-.191 


18 

0 

.026 

0 

.026 

20 

.090 

.018 

.090 

.018 

22 

.151 

.035 

.151 

.035 

24 

0 

0.028 

0 

0.028 

25 

-.128 

.020 

-.128 

.020 

26 

-.108 

.025 

-.108 

.025 

27 

-.103 

.014 

-.103 

.014 

28 

.175 

.020 

.175 

.020 

30 

.121 

.017 

.121 

.017 

32 

0.023 

0.024 

-0.241 

0.049 

34 

.058 

.020 

.211 

.026 

37 

.652 

.042 

.050 

.030 

40 

.552 

.045 

-.114 

.042 

43 

.579 

.035 

.148 

.040 

46 

.422 

.023 

.421 

.020 

49 

0.030 




52 

.908 



.039 

55 

.556 


-.156 

.027 

58 

.607 


-.159 

.035 

61 

.758 


.711 

.059 

64 

.438 


.198 

.073 

67 

-0.017 

0.018 

-0.245 


70 

.070 

.017 

.196 


73 

.461 

.022 

.100 


76 

.179 

.013 

-.085 


79 

.002 

.015 

.022 



11 



Table V. Wind-Tunnel Estimates for Lateral-Directional Stability Derivatives 


a, deg 

C Yp 

C b 

C b 

Crip 

0 

-1.150 


■ 

0.129 

5 

-1.200 


1 

.206 

10 

-1.150 


! 

.229 

15 

-.974 

-.103 

-.300 

.218 

20 

-.630 

-.160 

-.050 

.155 

25 

-.401 

-.160 

0 

.072 

30 

-0.230 


■ 

0.049 

35 

-.115 

■ 

1 


40 

-.372 

-.185 


1 

45 

-.343 

-.281 

.240 

.018 

50 

.057 

-.270 

.340 

.069 

55 

.343 

-.206 

.190 

.095 



-0.215 



65 

.286 

-.143 

(a) 

.329 


.143 

-.135 

(a) 


75 

.229 

-.143 

(a) 

.227 


.573 

-.149 

(a) 

.441 


“Data not available. 


12 










Figure 1. Radio-controlled X-29A drop model. 



Figure 2. Schematic representation of drop test procedure. 
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